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Abstract. We study the kinetics for the search of an immobile target by randomly 
moving searchers that detect it only upon encounter. The searchers perform 
intermittent random walks on a one-dimensional lattice. Each searcher can step on 
a nearest neighbor site with probability a, or go off lattice with probability 1 — a to 
move in a random direction until it lands back on the lattice at a fixed distance L away 
from the departure point. Considering a and L as optimization parameters, we seek 
to enhance the chances of successful detection by minimizing the probability Pn that 
the target remains undetected up to the maximal search time N. We show that even 
in this simple model a number of very efficient search strategies can lead to a decrease 
of P/v by orders of magnitude upon appropriate choices of a and L. We demonstrate 
that, in general, such optimal intermittent strategies are much more efficient than 
Brownian searches and are as efficient as search algorithms based on random walks 
with heavy-tailed Cauchy jump-length distributions. In addition, such intermittent 
strategies appear to be more advantageous than Levy-based ones in that they lead 
to more thorough exploration of visited regions in space and thus lend themselves to 
parallelization of the search processes. 
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1. Introduction 

Search processes are ubiquitous in nature: Predators search for prey, prey also hunt, 
molecules search for each other to recombine in order to produce required chemicals, 
proteins search for target sequences on DNAs. Human beings spend their lives searching 
for different things - better jobs, shelters, partners, lost keys; they also seek efficient 
search strategies to minimize the time to reach the desired target, or at least to enhance 
their chances of eventually finding it. 

The search for a desired target may depend on a variety of different conditions and 
may take place in different environments. Targets may be sparse, hidden, difficult to 
detect even when found. The targets may be mobile or immobile, they may try to avoid 
searchers, there may be one target or many. They may have a finite life-time and vanish 
before they are detected. Searchers, on the other hand, may search "blindly," detecting 
the target only upon encounter, or they may perceive distant targets and adjust their 
motion accordingly. They may have no memory of previously visited areas, or they may 
avoid such areas. The searchers may act individually or in swarms, optimizing their 
search efficiency by exchanging information. Finally, the "efficiency" of a search may 
be judged by a variety of measures, including the time to reach a target or targets, the 
number of encounters of searchers and targets per unit time, or the exploration range 
of space per unit time. In general, for each specific situation different search strategies 
may be appropriate. The quest for optimal strategies has motivated a great deal of work 



in the past few years [IlElElllElEllTlElElIIOliniliaiTailSlIiniini 
[221 [231 [22 [251 [261 [271 Eg. 




Earlier work tended to focus on deterministic search algorithms (see, e.g., Refs.[Tl[2l 
[31 [1] and references therein) specific to human activities such as, say, rescue operations 
or the search for natural resources. More recent studies have focused on random 
search strategies. In this context, it has become quite clear that strategies based 
on Levy flights (instantaneous) or Levy walks (occurring with a finite velocity) are 
according to all measures more advantageous than strategies based on conventional 
diffusive motion or random walks with steps to nearest neighbors only (Brownian 
search) O [6l [71 [H [9l [TOl [HI [12]. Levy searchers perform excursions whose lengths / 
are random variables with heavy-tailed distributions p{l) of the form 



where i? is a normalization constant. In particular, it was demonstrated in Ref. [29] 
that using the Cauchy distribution of jump lengths as given in equation ([1]) with /i = 1 
instead of a Gaussian distribution (Brownian search) allows for a faster cooling scheme, 
and hence for a considerable reduction of computer time in the search for a global 
minimum in nonconvex (multiple extrema) energy landscapes by simulated annealing. 
Aside from this, extensive data has been presented allegedly supporting the idea that 
many of species in the living world do indeed follow Levy-type random motion in their 
search [HI [6l [71 [H [9l [TOl [TTl [12]. This point, however, has been questioned recently in 
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Refs. [201 EI]- The main objections of Refs. [301 EI] have recently been re-examined in 
Ref. [12], where it was concluded that Levy-based strategies may still be consistent with 
experimental observations if one takes into account a highly non-homogeneous spatial 
distribution of prey. More sophisticated models with adaptive behavior in which the 
foragers use their cognitive skills to develop more efficient foraging strategies have also 
started to appear in the literature (see, e.g., Ref. p^]). 

Following the observation of trajectories of foraging animals such as lizards or fish 
or birds, in which active local search phases randomly alternate with relocation phases 
(see, e.g., Refs.|14[ [15]). another type of random search - an intermittent search - has 
been proposed. In this algorithm, a search is characterized by two distinct types of 
motion, a ballistic relocation stage during which the searcher is non-receptive to the 
target, and a relatively slow phase with random Brownian-type motion during which 
the target may be detected [161 Il3 UHl UHl EI] (see also Ref. [20j). In this approach 
one aims to minimize the time of first passage to the target from a given location by 
varying, for example, the relative durations of the active and relocation stages. 

In this paper we pursue the optimization of an intermittent strategy. A simpler 
version of our model was presented in Ref. [2T] , where we developed a search algorithm 
based on intermittent random walks that involve nearest neighbor steps and off-lattice 
relocations of a fixed length L. There, instead of minimizing the first passage time to 
the target, we sought to maximize the success of the search by minimizing the target 
non-detection probability over a fixed finite maximum search time. We note that the 
time derivative of this probability defines the distribution function of the first passage 
time to the detection event. Thus, contrary to previous work in which only the mean 
first passage time was optimized [161 [13 [HI [El [21] , our goal was to optimize the full 
distribution function. Our optimization parameter was the intermittency parameter a 
that determines whether the next step will be on- or off-lattice. It was shown, both 
analytically and numerically, that the probability of failure to detect the target over a 
finite search time can be made smaller by many orders of magnitude upon an appropriate 
choice of this parameter. We note that in Ref. [22] a different intermittent search 
algorithm was proposed in which the length of the relocation stage is a random variable 
with a heavy-tailed distribution in equation ([T]). In that work it was concluded that 
such a combined strategy is advantageous over intermittent searches with exponentially 
distributed [HI [13 [H] or fixed [21] relocation lengths since it allows a searcher to 
find the target more quickly in the critical case of rare targets, and since the search 
performance is much less dependent on adaptation to the target density. However, 
we argue that in some (albeit not all) situations, strategies involving Levy distributed 
relocations can not be optimal. This occurs when there is some maximal time that the 
search is allowed to run. Allowing the length of the relocation stage (and consequently, 
the time spent on each relocation tour) to have a heavy-tailed distribution would lead 
to some portion of the search process that would involve unnecessarily long relocations 
divorced from the targets, thus not contributing to the overall finite time effort. Thus, 
when appropriate, an optimal search algorithm should be based on relocation length 
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distributions that explicitly account for the fact that a search process is limited in time. 
It may well be that the optimal jump-length distribution should itself vary with time. 
Furthermore, trajectories of Levy walks or flights are "overstretched" in the sense that 
such walks explore space in a very irregular manner. The visited area consists of a 
patchy set of disconnected clusters, leaving large unexplored voids compared to the case 
of a Brownian search. Additionally, when many Levy searchers are involved, the fact 
that a Levy distribution does not have moments induces rapid mixing. This mixing 
might be advantageous if the detection probability is low (or the false alarm probability 
high), such that multiple rechecking of each site by other searchers is required in order 
to spot the target. Otherwise, this very efficient mixing might be a disadvantage since it 
does not favor the parallelization of the search process by dividing the searched area into 
subunits and assigning a separate domain to each searcher. In the "living" world animals 
are often constrained to their assigned territories, and even an occasional incursion into 
a neighbor's terrain while searching for prey may cause serious difficulties. 

In this paper we revisit the question of an optimal jump- length distribution 
underlying an efficient search algorithm. Focussing on the one- dimensional case 
(for which Levy-based search strategies are said to most dramatically outperform 
intermittent ones), we study the search kinetics of a "hidden" immobile target located at 
the origin of an infinite lattice by a concentration p of randomly moving searchers. The 
motion of the searchers is intermittent, consisting of two distinct, randomly alternating 
stages - ballistic, off-lattice relocations with finite velocity over a fixed distance L, and 
random walks between nearest-neighboring sites. In other words, we consider a search by 
intermittent random walkers with a jump-length distribution of the form [very different 
from that in equation ([T])] 

pil) = f m - 1) + + 1)] + m -L)+ 6{i + L)] , (2) 

such that the searchers step on nearest neighbors with probability a and, with 
probability 1 — a, perform long jumps over a distance L |2T]. The process evolves 
in discrete time n = 0, 1, 2, . . . , A^, where N is the maximal time the search process may 
run. This time may depend, for instance, on our patience or on experimental constraints. 
Note that the constraint of the maximal search time N is the crucial aspect of our work 
which makes our analysis very different from other models. Steps to nearest neighbors 
take unit time, while off-lattice relocations over a distance L require time T. The term 
"hidden" means that a searcher can not perceive the target when off-lattice. A searcher 
only detects the target when it arrives at the site on which the target is located. 

We pose the following question: Is it possible to choose a = aN and L = Ln, 
dependent on the maximal search time A^ but independent of the running time n, 
which optimizes the search efficiency and leads to a performance that is better than 
Levy-based strategies? In order to answer our question, our first goal is to calculate the 
probability P/v that the target remains undetected up to the maximal search time A^ and 
to determine its asymptotic behavior analytically in the large- A^ limit. Then, considering 
a = apf and L = L^r as optimization parameters, we seek to enhance the searchers' 
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chances of success by minimizing P^. We will demonstrate that, depending on whether 
we are at liberty to tune a (as in Ref. [2T]), or both a and L, different optimal strategies 
can be realized all of which can decrease the value of P/v by many orders of magnitude 
compared to a Brownian search. We also show that even the simple distribution ([2]) 
with optimal a = and L = yields very efficient search algorithms comparable to 
and in some cases better than Levy-based strategies. Moreover, in striking contrast to 
the latter, optimal intermittent walks lead to much denser exploration of space. These 
results support our claim that Levy-based searches are not the best algorithms when 
the search is limited in time. 

2. Model and basic equations 

Consider a one- dimensional regular lattice of unit spacing containing M sites labeled 
by s. The lattice is a circle, that is, we use periodic boundary conditions. At one of the 
lattice sites, say at the origin s = 0, we place an immobile target. Then we randomly 
place K searchers under the constraint that none is placed at the site of the hidden 
target. We focus on the behavior in the thermodynamic limit, K — * oo, M —>■ oo, with 
a finite mean density of searchers p = K/M. We note that the model under study can 
also be solved in the general case of finite K and M, but the calculations become more 
involved without adding significant new features to our conclusions. 

Next we allow the searchers to move according to the following rule. At each tick of 
the clock, n = 1,2,3, N, each searcher selects randomly between two possibilities: 
with probability a, it moves with equal likelihood to one of its nearest neighboring sites, 
and with probability (1 — a) it leaves the lattice and flies off-lattice with a given velocity 
V until it lands L sites away from the departure site. The distance L is fixed, but 
either direction of the flight is chosen at random with equal probabilities. The time it 
spends off- lattice during the flight is T = L/V. This parameter can take integer values, 
T = 1,2, . . . , L. Note that this condition defines the velocity V. There is no restriction 
on multiple searcher occupancy of the sites. Note as well that in the two "pure" cases, 
a = 1 and a = 0, the model reduces to standard random walks. Here we take both 
a and L to be independent of n but possibly dependent on the maximal search time 
N. An interesting situation with time-dependent a = a„ and L = will be discussed 
elsewhere [33]. We focus on perfect detection, that is, a searcher recognizes the target 
immediately upon first contact. The case of imperfect recognition can be solved using 
the same techniques and will also be discussed elsewhere [33] . 

The probability that the target has not been detected by step N is related in a 
simple way to Sn-, the number of distinct sites visited up to that time, as 

Pn = exp {-pSn) (3) 

(see, e.g., Ref. [35]). This result, which holds for independent searchers and was 
explicitly shown to be valid for our model in plj, is a crucial equation since we will 
arrive at results for P/v via calculations of S'at. A larger leads to a smaller probability 
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that the target remains undetected and thus to a better search algorithm. In general, 
Sn ~ AN'' with < 7 < 1. One thus expects that larger 7 leads to a more efficient 
search and explains why, intuitively, it was believed that the most efficient search 
algorithm should be based on Levy walks with a broad distribution of jump lengths, for 
which Sn grows more rapidly than in the case of simple Brownian motion. As an aside, 
however, we note that even for standard Brownian motion S^- ~ AN/ In(A^) in two 
dimensions and ~ AN in three dimensions; consequently, choosing Levy walks as a 
search mechanism will not lead to any significant gain compared to a Brownian search 
in these higher dimensions (except perhaps through the prefactor A). In one dimension, 
however, there are significant differences in the growth rates of Sn between Levy and 
Brownian motions, and our task is to explore the place of intermittent random walks in 
this panorama. 



3. Expected number of distinct visited sites 

The expected number of distinct sites visited can be calculated once we determine the 
probability P{s\sQ;n) that a given searcher, starting its intermittent random walk at 
site So at time moment n = 0, appears at site s at time moment n. More specifically, 
the generating function S{z) of Sjy, defined as 

00 

S{z) = Y: Snz"", (4) 

N=0 

and the lattice Green function (or site occupation generating function) 

00 

P{s\so;z) = Y.P{s\so;n)z'' (5) 

n=0 

of the intermittent random walk are related to each other through [21] 

^^'^-l-z P(0|0;^) • 

Hence, given P{s\0;z), we obtain S{z) by virtue of equation (Q. We then determine 
the iV-dependence of Sn by inverting the discrete Laplace transform in equation (jl]). 
The probability P(s|so; n) obeys the recurrence relation 

<y 

P{s\so; n) = - [Pis - l|so; n - 1) + P(s + l|so; n - 1)] + 
1 - a 

+ [Pis - L\so; n-T) + Pis + L\so; n - T)] , (7) 

which explicitly takes into account that jumps between nearest-neighboring sites proceed 
in unit time, while long-range jumps over distance L require an integer time T. 
Equation ([7]) thus defines a non-Markovian process with a memory. Note also that 
since the intermittent random walk defined by equation ([7]) is homogeneous so that 
P(s|so; n) = Pis — so\0; n), without loss of generality we henceforth set sq = 0. 
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Multiplying both sides of equation ([7]) by z'" and performing the summation, we 
find that P{s\0] z) obeys 

1 cos (ksjdk 

pmz) = - , \/ (8) 



vr JO 1 — az cos{k) — (1 — a)z'^ cos{kL) 
and consequently the generating function of the expected number of distinct visited sites 
is given by 

' (9) 



z) {1 — az — {1 — a)z'^) 



dk 



1 — azcos{k) — (1 — a)z'^ cos{kL) 

Before we proceed further, the following remarks are in order. Note that P{s\0]z) in 
equation ([HD, and consequently S{z) in equation ([H]), can be calculated explicitly in the 
two "pure" random walk cases, a = 1 and a = 0. In these two limits one finds for 
the large- asymptotic behavior 

/Q]\f \ 1/2 / 1 \ 

M« = i)=( — ) (10) 



TT 



and 



Ma = o)^[-) +o[j^_). (11) 

The result in equation ffTOl) is well known (see, for example, Ref. |35j). 

• Sublinear growth of Sn with time signifies that each site visited by such a walk 
is most probably visited many times. This oversampling is precisely the reason why 
searching a target in a one-dimensional system by a Brownian search is not very 
efficient, since the walker wastes a great deal of time revisiting sites that do not 
contain the target. That is why, in fact, recourse has been made to Levy-based 
searches, since they reduce oversampling and lead to stronger growth of Sn with 
N. 

• The result in equation (fTTj) is the same as that in equation (fTOjl with the replacement 

N/T and hence does not represent a good search strategy either - in fact, it 
is worse. When relocations over a distance L take unit time as do nearest-neighbor 
steps, the results in equations ffTOjl and (|TT|) coincide, as they should. 

• Note as well that in one dimension, equation ([3]) with S^- in equation (fTOl) defines 
the asymptotically exact behavior of the non-detection probability of a target which 
dijfuses in the presence of a concentration of diffusive searchers [36]. Thus the 
asymptotic behavior of P/v is independent of the target diffusion coefficient (see 
also Refs. ^ and [23l El ESI EEl [27] for more details). 

• The result in equation ([3]) can be generalized to the case of imperfect recognition of 
the target, that is, when target recognition upon encounter occurs with probability 
p < 1 [34]. In one-dimensional systems the leading asymptotic behavior of P/v is 
independent of p provided that p > 0, and thus is also described by equation ([3]). 
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We seek a large- iV expansion of Sn, in which (to arrive at the correct optimization) 
it is essential to retain not only the leading divergent contribution as oo but 

also, if present, a constant A^-independent correction term. Turning to the limit 
2; — i> 1~ (A^ ^ 1) and inverting equation ([9]) we find, after some rather tedious but 
straightforward calculations, that Sn obeys 

SN = fiN'^' + f2 + 0(-^), (12) 



where the coefficient /i is given by 

1/2 

/i = ( -— ^ 1 , a > 0, (13) 




a = 0. (14) 



The parameter 

a 



, (15) 
I — a 

is an important physical parameter which defines a characteristic time for a "continuous 
tour of diffusion," that is, the typical time spent by a searcher on the substrate between 
two consecutive off-lattice flights. Note that the leading term in equation (fT2ll grows 
as A^^/^, which means that the leading behavior is that of a random walk, albeit 
intermittent, unless there is an additional dependence of /i on A^ via the optimization 
of equation ( |T2l) with respect to a and L. Note also that the coefficients /i, /2, • • • are 
discontinuous functions of a, and that a = is a singular point since it is not possible for 
a random walker that skips over L sites at each step to ever visit all sites, cf. Ref. |32j . 
This discontinuity should be viewed with appropriate caution since for any fixed finite 
A^, S]sf is a smooth function of a; the discontinuity arises because we are describing an 
asymptotic behavior that is strictly valid only in the N ^ 00 limit. Additional details 
and explanations can be found in Refs. [21] and [32] . 

Returning to the asymptotic A^ — »• cxd limit, for a = the constant term /2 = 0, 
while for < a < 1 (note that this double-sided inequality is strict) it is determined by 

2(a + (l-a)L2) fU.. du 



2{a + {V-a)L-) 
h — / — / 



n\/a{l-a) -^o sh.{u) ^ 1 - rsli^ {u 
sh^{2Lu) L 



2sh\Lu) + Tsh\u) T + L"^ 



where 



cth(M) , (16) 



Um = -\n\ + \{ -1 . (17) 




The integral in equation f|T6|) can not be performed in closed form, but its important 
contribution to the problem can be estimated. Anticipating that effective search 
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strategies take place when 2UmL ^ 1 (see below), we find that in this limit the leading 
behavior of /2 is given by 

h-~l{r'^' + ^^9iL), 0<a<l, (18) 

where g{L) is a slowly- varying function of L, 

Equations (ll2p - f[T^ constitute our main general result and provide the basis for the 
design of an optimal strategy through the choice of a and L. Below we discuss such 
a design and show that, indeed, the optimal search strategies fulfill the assumption 
2UmL > 1. 



3.1. Optimization: Tuning a at fixed L 

To highlight the optimization procedure, we start with the case studied analytically and 
numerically in [2T], shown again here for completeness, namely, we tune a but hold the 
relocation length L fixed. 

Note that Sn defined by equations (|T2l) - (fT9l) is a non-monotonic function of the 
characteristic diffusion time r. Differentiating with respect to r (discarding a weak 
dependence of g{L) on r), we find that the maximum of with respect to r is given 
implicitly as the solution of the equation 

^ ' ^ 'l- — U(L), (20) 



dr \ 71 J 7rr^/2 



where 



(21) 



dfi^l ( 1 {t + I^\ 

dr 2 l^[(r + T)(r + L2)]i/2 (r + T)3/2 j " 

The left-hand-side of equation fl20p diverges when — ^> oo, which indicates that in this 
case with fixed L the optimal time r of continuous tours of diffusion should tend to zero. 
We find that to leading order in the optimal r = Topt and hence, the optimal value 
aopt of the intermittency parameter a, are given by 

L2/3ln2/3(L) ^ ^ 

aopt ~ Topt ~ T ^^^j^y/s ■ (22) 

This is consistent with the condition 2UmL S> 1 since UmL ~ L\n{l/a) ~ Lln(A^) ^ 1. 
The symbol ~ here and henceforth signifies the exact behavior to leading order in N. 
The expected number of distinct sites visited by an intermittent random walk with an 
optimal a and fixed L then is 

S.-^C-^r, (23) 



TT 



i.e., it differs by a factor of L/y/T from the corresponding result for a standard nearest- 
neighbor random walk with a = 1 (Brownian search), equation fllOp . 
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The essential result of this subsection is an enhancement by a factor L/^/T of 
the expected number of distinct sites visited by an intermittent random walk with the 
distribution in equation ([2]) and an appropriate A^-dependent choice of the intermittency 
parameter compared to the outcome of an ordinary random walk. Note that this effect 
appears in an exponent in the non-detection probability, so it can become dramatically 
apparent. For example, for L = 5, T = 1 and = 10^, with a density of searchers 
as low as p = 0.01, the non-detection probability for a Brownian search is P/v ~ 0.2 
while that of the optimal intermittent search {aopt ~ 0.07 for these parameters) we 
find P/v ~ 0.0003, a reduction of three orders of magnitude! Note finally that in one 
dimension for fixed L the optimal strategy involves a progressively smaller fraction of 
nearest-neighbor steps as is increased. 



3.2. Optimization: Tuning a and L for T = 1 

Next we consider both a and L in equation ([2]) to be tunable, but we fix T = 1, that 
is, relocation to a nearest neighbor and to a neighbor a distance L away both take one 
unit of time. This causes the relocation velocity V = L/T to become dependent on N 
via L. Setting T = 1 in equations f|T2l) - f|T9l) . we have 



S„ = (a + (l-a)L^) (-) --((—) +(— ) L^)9(L). (24) 

Note that the first term on the right hand side of equation fl24l) grows with L and thus 
favors high values of L, but that the second term is negative and contains a higher power 
of L. This "competition" suggests that there exists an optimal A^-dependent value of L 
which leads to a maximum in the number of distinct sites visited. 

To determine the optimal value of L, we differentiate Sn with respect to L. Again 
discarding a logarithmically weak dependence of g{L) on L and anticipating that the 
optimal value of L is large, such that L ^ (ct/(l — «))^^^ (to be checked later for 
consistency), we find that the optimal value L = Lopt obeys 

L^^^MLopt) - [-^) , (25) 

and hence. 



Substituting this result into equation (12^ then leads to the expected number of distinct 
sites visited in an intermittent random walk with an optimal length Lopt of the relocation 
stage, 

S^^2{a{l~a)f \ \ (27) 

ln(7ra;A'/2) 

Note the parabolic form of the prefactor as a function of a. The prefactor vanishes in 
the pure limits a = and a = 1, indicating that optimization in these pure cases is not 
possible. This is a reflection of the fact that the leading behavior of in these limits 



Efficient search by optimized intermittent random walks 



11 



is determined by terms proportional to N^^'^ and constant, independent of terms are 
absent. 

Further optimizing the prefactor in equation ( 1271) with respect to a, we find that 
«opi — i> 1/2 as — cx), and hence 

^ ^^^^ 

This result is consistent with the earlier assumptions 2UmL 3> 1 and L ^ [a(l — a)]^^^ 
since here Um and a /{I — a) are constants while Lopt diverges as — oo. 

The result (128|) shows that when T = 1, optimization of the intermittent search 
with respect to both a and L leads to an additional factor N^^"^ / In(A^) in Sn via 
the coefficient fi in equation f|T2|) . which results in a much stronger dependence of 
the expected number of distinct sites visited on the maximal time A^. In fact, we 
have obtained a behavior close to that of a two-dimensional Brownian motion, which 
signifies that optimal intermittent random walks lead to only marginal oversampling. 
The optimal strategy here consists of taking aopt = 1/2 {jopt = T = 1), and relocation 
length Lopt as given in equation fl26|) . Note that a similar result, i.e., that the maximum 
Sn is attained when the time spent on relocations is equal to the time spent in the 
diffusive stage, has been obtained for a model describing a stochastic search of a target 
site on a DNA by a protein [3HI HH] . 

At this point one might be tempted to say that exactly the same temporal behavior 
of Sm as in equation fl28l) can be found without resorting to any optimization procedure 
but by merely taking a Levy walk with /x = 1 in equation ([T]) (Cauchy distribution). 
Indeed, in this case one obtains (see equation (2.20) in [39] ) 

where the normalization constant B in equation ([T]) with fi = 1 has been set to S/tt^. 
Remarkably, comparing the prefactors in equations fl28l) and fl29l) one notes that the 
search based on the intermittent strategy with fine tuning of the optimization parameters 
outperforms the one based on a Cauchy distribution due to a numerical factor which is 
more than six times larger in the intermittent walk. 

Moreover, we emphasize that these algorithms are very different in their quality of 
exploration of space. Consider, for example, the "density of visited sites" defined by 

= it- <'°' 

where Mn is the expected maximal displacement in, say, the positive direction so that 
2Mjv is a measure of the range of the walk. The parameter fijy is thus a measure of how 
many sites have been visited within the range of the walk. By definition, < f2Ar < 1. 
For a Levy walk with a Cauchy distribution of the relocation length the distribution of 
the maximal displacement is well-known [lO]; the tail of this distribution exactly follows 
the behavior of the parent Cauchy variables and hence, ]\^^°-^'^'^y jg infinite. This implies 
that the density of visited sites vanishes, fl^ = 0, and thus the exploration quality is 
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very poor. The expected maximal displacement of the intermittent random walk can 
be found from the general result of Ref . [H] , 



where 



1 /2 

2M^ = (a + (1 - ^V^) + In + O (-^) , (31) 



2 dk ^ /2(1 -acos(fc) - (1 -a)cos(fcL))\ , , 

^^ = ^70 P^H (a + (l-a)L2)P )■ (^'^ 

We thus find that fi^v ^ 1/2 as ^ oo, which signifies that intermittent random 
walks have a very good exploration quality in that they visit half of the sites within 
their range. This is not a result expected a priori since we are dealing with random 
walks that involve steps not only to nearest neighbor sites but also to distant sites. We 
note that the exploration quality may be further enhanced by optimizing both Sn and 



3.3. Optimization: Tuning a and L for fixed V 

We finally turn to the most difficult case, when relocation over distance L proceeds with 
a finite fixed velocity V. Note that this differs from the previous case, where T = L/V 
was fixed. We take note of two points: 

• We expect that aopt ~ 1 since fiights over distance L are less favorable now that 
each relocation costs time T = L/V during which no new sites are visited. If it 
turns out that the optimal relocation distance again grows with (as it does, see 
below), then the time T grows with as well. This in turn implies that it might 
become more advantageous to remain on the lattice, which means that the optimal 
r might be larger than in the previous case. 

• On the other hand, the expression in equation f|T8|) is only valid when 2UmL ^ 1. If 
dopt ~ 1, then Um ~ ^(1 — «)/« = l/y/r <^ 1. Consequently, approximation (ITSjl 
will be valid if the optimal characteristic time r is not too large, i.e., r -C L^. 

Differentiating equation ( fT2|) with respect to r (again discarding the logarithmically 
slow variation of g{L) with r), we have 

/ 1 VTTl? \ /8ivy/2 _ 2 / 1 

[^iT + L/V)iT + L^) (r + L/Vf/^) U J ~ r^^ 
Since equation (ITH]) is only valid when t <ti L^, equation (15^ can be simplified to yield 
L f8N\'/^ 2 L2 

Now we proceed as follows. We first assume that r ^ L/V, find Topt and L^pt, and 
ascertain whether the assumption is valid. Then we will follow with a more accurate 
calculation. 




Efficient search by optimized intermittent random walks 



13 



If r < L/V, then 

°^*"y(2vriV)V3- ^^^^ 
Substituting this expression into equation (fT2l) . we have 

V TT / TT 

Differentiating the latter expression with respect to L, we find that Lopt is defined 
imphcitly by 

L%',g'/'{Lopt) = ^{27iNf\ (37) 
which imphes that 

Hence, Xop^ is not much smaller than Lopt/V . It is only smaller by a numerical factor 
and scales as Lopt ~ N^^"^ / hi{N) (note that nonetheless Topt -C -^^^pj. Note also that 
Topt — >■ cxD as — > oo, which implies that Uopt 1 and Um ~ (1 — a)'^^"^ ~ r^/^. 

We next try to search for optimal L and r from equation fl34|) supposing that 
Um ~ T^''^ and r = CL/V ^ where C is a constant to be determined. We note that in 
this case 

g{L) = ln(L) - In ((1 + r)Vl - «) + 0.126 + O Q 
^hn{L)-^ln{C/V) 

^lln(L). (39) 

Hence, equation (134|1 becomes 

\/3/2L /8iV\^/2 2 L2 1 



In (L) . (40) 



(41) 



(1 + C)3/2L3/2 V TT ; TT T^/S 2 

From equation (HOl) we find that the optimal value of r obeys 
_ (1 + C) In'/^ (L) 
"""^^ " (27riV)^/^ ' 

and ^Af, optimized with respect to r, then follows: 

Differentiating this equation with respect to L, we find that the optimal flight length L 

obeys 

.3.3/2 



Lopt 



In (L„p,) = (^^ j (87riV)^/\ (43) 
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so that 



(IT (s^nY 

In ((f) (SttAT)^/^ 



Lastly, we obtain C from the definition Topt = CLopt/V. This gives 
n ^ -^opt In^^^ (-^) _ ^^opi 

^' + ^^1^ (SvriV)^/^ "^"^ ^^^^ 
which is solved by 

C = I (46) 

Thus the optimal strategy is realized when we choose Lopt as in equation (l44l) . and 
r = (3/4)T, i.e., for the optimal strategy the characteristic time of a tour of diffusion 
between two consecutive long jumps is three-fourths of the time a searcher spends on 
jumps over a distance L. 

Finally, we combine these results to find that the expected number of distinct sites 
visited optimized with respect to both a and L at fixed V is 

Note that this more intricate optimization procedure results in a faster growth law of 
Sn with than in the pure random walk cases. Consequently, the search efficiency has 
again been enhanced by orders of magnitude, albeit not as much as in the fixed-T case 
of the previous subsection. 



4. Conclusions 



We have considered the design of an optimal search strategy of a hidden target by 
a given density p of random walkers who have a limited maximal time to find 
the target. The analysis is restricted to one-dimensional systems. Our measure for 
the quality of a strategy is the minimization of the probability P/v that the target is 
undetected within the given maximal search time. In particular, we consider strategies 
that consist of a combination of nearest neighbor walks and jumps of fixed length L, both 
involving steps in either direction with equal probability. The motion of each searcher 
is thus intermittent. The probability that the target is undetected is just the survival 
probability for the target and is related to the distinct number of sites Sn visited by a 
walker up to time by the well-known relation = exp(— p^at). Our goal has thus 
been to maximize Sn- We stress that the time derivative of P/v defines the distribution 
function of the first passage time to the detection event. This means that in contrast 
to previous work |16l [TTl [18], [191 [21], our aim has been to optimize the full distribution 
function and not only its first moment. 

Our model has three parameters: a, the probability that the next step of the 
walker is a nearest neighbor step; L, the length of a long step; and T = L/V, where 
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T is the time it takes to cover a long step and V is the velocity of a long step. The 
parameters a and L are optimized as a function of the maximal time under different 
constraints. We compare our results for Sn with those of a nearest neighbor random 
walk, Sn ~ (8A^/vr)^/^, and of a Levy walk with a Cauchy distribution of step lengths, 
Sn ~ {3/27r^)N/\n{N). 

If L and T are fixed and only a is picked for optimal strategy, the best choice is to 
take it to be very small, a ~ A^~^/^ [21]. Most of the random motion then consists of 
steps of length L. The resulting Sn is larger than that obtained with a nearest neighbor 
random walk by the numerical coefficient L/^/T, Sn ~ (L/T^/^)(8A^/7r)-'^/^, so that the 
non-detection probability, which involves this factor in the exponent, can be decreased 
by orders of magnitude even for a low density of searchers. 

If we optimize both a and L, the best choice of these parameters depends on the 
constraint we place on the third parameter. It the time for a long step is the same 
as the time for a nearest neighbor step, then the optimal choices are a = 1/2 and 
L ~ N^/"^ /ln{N). Thus, short and long steps should occur with equal likelihood, and 
the optimal distance covered by the long steps grows (slowly) with increasing observation 
time. In this case the distinct number of sites visited is even larger (again by a numerical 
factor) than that of a Levy walk, Sn ~ N/ ln('7rA^/4). These N dependences are the same 
for any fixed value of T, although the specific numerical coefficients depend on this value. 
We find an important difference between the Levy walk and our optimized intermittent 
walk in the coverage of space, which may be an important feature if one wishes to 
parallelize the searches of different walkers. The density of visited sites vanishes in the 
Levy case, while that of the intermittent walk approaches a constant with increasing A^. 

Finally, if we again optimize both a and L but now keeping the velocity V of the 
long steps fixed (which means that the time for a long step grows with the length of the 
step), we find the optimal choice a to be close to unity, a ~ l-4Vln{N)/3N'^/^, and 
the optimal length step to grow with as L ~ N^^"^ /ln{N). Now the walkers rarely 
jump over long distances, but the time spent diffusing and the time it takes to make a 
jump of length L both grow with increasing observation time N , so that r = 3T/4. The 
distinct number of sites visited here shows a growth intermediate between the other two 
intermittent strategies, Sn ~ V^^/^A^^/^/ln^/^(A^). 

We stress that in all the strategies considered here, we have implemented a 
maximum observation time as part of the optimization process, and we have used a 
particular measure of the quality of a strategy, namely, that the probability of survival 
of the target be minimized within this time. We have shown that in one dimension 
even a simple intermittent step distribution consisting of nearest neighbor steps and 
long steps of fixed A^-dependent length can be optimized so as to yield far better search 
outcomes than a nearest neighbor random walk and even a Levy walk with a Cauchy 
distribution of step lengths. The model can be generalized and further optimized in a 
number of ways, some of which we have noted in the course of our analysis. 
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